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1.0  INTRODUCTION 


1.1  PROBLEM  DEFINITION 

Turbulent  near  wakes  behind  blunt  bases  in  supersonic  external  flow  occur  in  a 
variety  of  flows  of  practical  interest.  The  near-wake  flow  process  (Fig.  1 )  is  characterized 
by  a  strong  interaction  between  the  inviscid  external  flow  and  the  turbulent  portion  of 
the  flow;  therefore,  computation  of  the  entire  flow  field  requires  a  coupled  solution  for 
the  viscous  and  inviscid  regions. 

Several  features  of  the  near-wake  flow  field  are  shown  in  Fig.  1 .  First,  a 
recirculation  region  exists  just  downstream  of  the  base.  With  no  bleed  flow  from  the 
base,  the  velocity  on  the  centerline  is  zero  at  the  base  plane.  When  flow  is  bled  into  the 
base,  the  onset  of  recirculation  is  delayed  until  the  bleed  flow  has  been  entrained  into 
the  turbulent  shear  layer.  If  the  bleed  flow  is  sufficiently  large,  the  recirculation  region 
can  be  eliminated.  As  shown  in  Fig.  1,  the  pressure  just  downstream  of  the  base  is 
considerably  lower  than  the  pressure  on  the  body  surface  upstream  of  the  base.  The 
centerline  static  pressure  rises  slightly  with  distance  until  the  bleed  flow  is  completely 
entrained,  then  tends  to  decrease  slightly  with  distance  as  the  recirculation  velocities 
become  larger.  At  some  distance  downstream  of  the  base,  the  shear  layer  merges  with  the 
axis  and  the  flow  begins  to  recompress.  The  recompression  region  is  characterized  by 
large  pressure  gradients  and  by  acceleration  of  the  flow  along  the  centerline.  Downstream 
of  the  rear  stagnation  point  (RSP)  the  velocities  throughout  the  viscous  layer  are  positive. 

If  the  flow  is  planar,  and  if  the  external  flow  upstream  of  the  base  is  effectively 
infinite  and  parallel  to  the  centerline,  then  the  centerline  static  pressure  in  the 
recompression  region  monotonically  rises  to  the  free-stream  static  pressure.  However,  in 
an  axisymmetric  flow,  the  pressure  in  the  fecompression  region  "overshoots"  the 
free-stream  pressure  (Fig.  1).  Downstream  of  the  recompression  region  the  flow  relaxes 
toward  a  self  preserving  turbulent  far-wake  structure,  in  which  static  pressure  gradients 
are  negligible. 

Experience  has  shown  that  the  base  pressure  and  the  near-wake  flow  structure  are 
dependent  upon  (1)  the  thickness  and  shape  of  the  boundary-layer  velocity  profile  just 
upstream  of  the  base,  (2)  the  inviscid  external  flow  properties  (Mach  number,  flow  angle, 
etc.)  at  the  base  plane,  and  (3)  the  gas  properties  and  flow  rate  of  the  base  bleed.  The 
near-wake  flow  structure  also  depends  on  the  lateral  extent  of  the  supersonic  external 
flow.  If  the  external  flow  is  bounded  by  a  solid  wall  or  by  a  free  boundary,  then  waves 
which  are  generated  by  the  expansion  process  at  the  base  plane  will  be  reflected  from  the 
boundary.  These  reflected  waves  can  influence  the  near-wake  structure  if  they  intersect 
the  viscous  wake  upstream  of  the  point  where  the  flow  on  the  axis  has  reaccelerated  to 
supersonic  speed. 
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Figure  1.  Features  of  the  near-wake  flow  field. 
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Several  examples  of  supersonic  near-wake  flows  which  are  of  interest  at  the  AEDC 
are  shown  in  Fig.  2.  The  classic  turbulent  near-wake  flow  is  that  behind  bodies  in 
supersonic  flight  (Fig.  2a).  Truncated  plug  exhaust  nozzles  (Fig.  2b)  have  been  used  for 
both  turbine  engines  and  rockets.  An  accurate  prediction  of  the  plug  base  pressure  is 
required  for  the  analytical  assessment  of  the  thrust  performance  of  such  nozzles.  The 
supersonic  external  flow  is  limited  in  lateral  extent,  and  wave  reflections  from  the  free 
boundary  must  be  taken  into  account  in  the  prediction  of  the  near-wake  structure. 
Truncated  plug  nozzles  are  often  operated  with  base  bleed. 

The  annular  ejector  (Fig.  2c)  is  often  used  in  place  of  the  traditional  ejector,  which 
has  a  centrally  located  supersonic  driving  nozzle.  The  supersonic  external  flow  in  the 
annular  ejector  is  bounded  by  a  solid  wall,  and  wave  reflections  from  the  boundary 
strongly  interact  with  the  near-wake  structure. 

A  highly  complex  near-wake  flow  structure  is  involved  with  the  external  burning 
propulsion  concept  (Fig.  2d)  which  was  first  proposed  by  Strahle  (Ref.  1).  According  to 
Strahle,  waves  are  generated  by  the  mixing  and  combustion  process;  these  waves  interact 
with  the  near  wake  and  cause  the  base  pressure  to  be  much  higher  than  that  with  no  fuel 
injection.  Indeed,  Strahle  suggests  that  the  base  pressure  can  be  increased  to  a  level 
greater  than  free-stream  static  pressure,  perhaps  to  a  level  where  net  thrust  can  be 
obtained  with  a  slender  forebody.  This  external  burning  propulsion  concept  has  been  the 
subject  of  several  experimental  investigations  (e.g..  Ref.  2).  It  is  clear  that  an  adequate 
analysis  of  this  complex  flow  must  be  based  on  a  physically  perceptive  analysis  of  the 
simpler  flow  shown  in  Fig.  2a. 

In  view  of  the  numerous  practical  flows  in  which  supersonic  near-wake  mechanisms 
occur,  it  was  decided  that  a  generalized  analysis  for  such  flows  was  needed  at  the  AFDC. 
The  generalized  analysis  should  include  the  following  factors:  (1)  arbitrary  turbulent 
initial  boundary  layers;  (2)  arbitrary  inviscid  external  flow  profiles  at  the  base  plane;  (3) 
external  flows  which  are  bounded  by  a  solid  wall  or  by  a  constant  pressure  boundary;  (4) 
arbitrary  base  bleed,  including  gases  which  can  react  chemically  with  the  external  flow; 
and  (5)  planar  or  axisymmetric  geometry.  The  analysis  should  provide  reasonably 
accurate  predictions  not  only  of  the  gross  features  of  the  flow,  such  as  the  base  pressure, 
but  also  of  the  detailed  flow  structure  in  the  viscous  wake  and  in  the  adjacent  inviscid 
flow  field. 
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1.2  REVIEW  OF  EARLIER  ANALYTICAL  WORK 

The  analysis  of  turbulent  near- wake  flows  has  received  much  attention  during  the 
past  twenty-five  years.  Most  of  the  work  has  been  based  either  on  the  viscous  interaction 
model  proposed  by  Crocco  and  Lees  (Ref.  3)  or  on  the  component  model  first  developed 
by  Korst  for  turbulent  flows  (Refs.  4  and  5).  Most  investigators  have  concentrated  on 
planar  flows,  which  are  inherently  simpler  than  axisymmetric  flows. 

Crocco-Lees  Model  -  This  integral  method  is  based  on  solving  the  integral  form  of 
the  boundary-layer  equations  for  the  viscous  portion  of  the  flow,  including  the  reversed 
flow  region.  The  planar  external  flow  field  is  coupled  to  the  viscous  region  by  the  use  of 
the  Prandtl-Meyer  equations,  which  relate  the  local  static  pressure  to  the  flow  angle  at 
the  edge  of  the  viscous  layer.  The  viscous  layer  growth  rate  is  computed  with  an 
empirical  entrainment  coefficient.  A  major  problem  in  the  original  Crocco-Lees  analysis  is 
establishing  the  relationship  of  the  integral  parameters  of  the  viscous  flow  profiles  to  the 
other  variables  of  the  problem. 

A  major  feature  of  the  Crocco-Lees  solution  is  the  appearance  of  a  saddle-point 
singularity  (the  critical  point)  in  the  recompression  region.  This  singularity  is  interpreted 
as  being  analogous  to  choking  in  the  throat  of  a  nozzle,  i.e.,  the  wake  flow  upstream  of 
the  singularity  is  subcritical  and  can  transmit  disturbances  in  the  upstream  direction.  The 
wake  flow  downstream  of  the  singularity  is  supercritical,  i.e.,  supersonic  on  the  average. 
The  correct  near-wake  solution  is  taken  as  that  which  will  pass  through  the  singularity. 

Although  applicable  to  a  broad  class  of  strong  viscous  interaction  problems,  the 
original  Crocco-Lees  analysis  yields  rather  poor  results  for  planar  turbulent  near  wakes. 
The  first  successful  extension  of  the  Crocco-Lees  framework  to  the  planar  turbulent 
near-wake  problem  was  made  by  Alber  (Refs.  6  and  7).  Alber  obtained  the  integral 
parameters  of  the  viscous  flow  profiles  by  using,  in  a  transformed  incompressible  flow, 
the  Stewartson  solutions  of  the  Falkner-Skan  equations.  The  empirical  mixing  coefficient 
of  the  original  Crocco-Lees  analysis  was  replaced  by  the  solution  of  the  first  moment  of 
momentum  equation  (the  mechanical  energy  equation);  the  dissipation  integral  in  this 
equation  was  evaluated  by  the  use  of  an  empirical  eddy  viscosity  model.  Prandtl-Meyer 
coupling  with  the  external  flow  was  retained, 

It  is  only  recently  that  an  attempt  has  been  made  to  extend  the  Crocco-Lees 
framework  to  axisymmetric  turbulent  wake  flows.  Strahle  and  his  coworkers  (Refs.  8  and 
9)  have  been  developing  an  axisymmetric  version  of  Alber 's  analysis,  with  the  objective  of 
extending-  it  to  the  external  burning  populsion  problem.  The  external  flow  is  assumed  to 
be  isentropic  and  is  treated  with  an  approximate  method  of  characteristics  technique. 
(Prandtl-Meyer  coupling  is  not  applicable  to  axisymmetric  flow).  Only  flows  without  base 
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bleed  are  considered.  The  solution  for  only  one  axisymmetric  flow  field  has  been 
presented  (Ref.  9). 

Korst  model  -  The  Korst  component  model  (Ref.  5)  is  much  simpler  than  the 
Crocco-Lees  model  and  can  be  considered  a  first-order  solution.  The  component  processes 
involved  in  the  near-wake  flow  field  (inviscid  external  flow,  turbulent  mixing, 
recompression)  are  considered  separately.  The  results  from  the  component  analyses  are 
then  combined  to  yield  a  solution  for  the  base  pressure  and  for  the  other  global 
properties  of  the  recirculating  flow  region.  Three  major  assumptions  are  involved  in  the 
basic  Korst  analysis:  (1)  the  recirculation  region  is  quasi-stagnant,  (2)  the  flow  along  the 
"discriminating"  streamline  (that  which  separates  the  flow  which  passes  through 
recompression  from  that  which  is  recirculated  upstream)  isentropically  recompresses,  and 
(3)  the  flow  along  the  discriminating  streamline  stagnates  at  the  peak  of  the 
recompression  pressure  distribution.  Although  experimental  evidence  is  available  to  show 
that  the  first  and  third  of  these  assumptions  are  quite  incorrect,  the  basic  Korst  analysis 
(especially  the  simplest  version  in  which  initial  boundary-layer  effects  are  neglected) 
yields  surprisingly  good  predictions  of  the  base  pressure  for  planar  and  quasi-planar  flows. 
The  analysis  has  been  extended  to  include  base  bleed,  asymmetrical  two-stream  flows, 
nonadiabatic  flows,  and  flows  with  chemical  reactions  in  the  near  wake  (Ref.  10). 

Numerous  attempts  have  been  made  to  improve  the  component  analyses  in  the 
Korst  model.  Nash  (Ref.  11)  introduced  an  empirical  factor  to  account  for  the  fact  that 
the  discriminating  streamline  does  not  stagnate  at  the  peak  of  the  pressure  distribution. 
Other  investigators  (e.g..  Lamb  and  Hood,  Ref.  12)  have  developed  control  volume 
analyses  for  the  recompression  region. 

The  basic  assumptions  in  the  basic  Korst  analysis  are  not  very  realistic.  Because  of 
compensating  effects,  the  simplest  Korst  analysis  for  negligible  initial  boundary  layers 
yields  fairly  satisfactory  base  pressure  predictions  for  experimental  flows  with  initial 
boundary  layers  of  moderate  thickness  (5  less  than  the  base  height).  Attempts  to 
properly  include  the  influence  of  the  initial  boundary  layers  in  the  Korst  analysis,  using 
techniques  such  as  the  apparent  mixing  layer  origin  shift  concept  of  Hill  and  Page  (Ref. 
13),  have  not  been  very  successful.  However,  these  methods  work  reasonably  well  when 
applied  to  a  mixing  process  which  does  not  involve  a  large  expansion  at  the  base  plane. 
The  discrepancy  for  near-wake  flows  is  related  to  the  fact  that  the  origin  shift  concept 
includes  the  change  in  the  boundary-layer  mean  flow  profile  when  it  passes  through  the 
corner  expansion  process,  but  does  not  take  into  consideration  that  the  boundary-layer 
turbulence  decays  rapidly  when  it  passes  through  the  expansion.  This  rapid  decay  of  the 
turbulent  fluctuations  in  the  "remnant"  of  the  initial  boundary  layer  was  observed 
experimentally  by  Lewis  and  Behrens  (Ref.  14)  and  Page  and  Sernas  (Ref.  15).  These 
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observations  lead  to  the  conclusion  that  the  boundary-layer  flow  after  the  corner 
expansion  should  not  be  taken  as  the  initial  condition  for  the  subsequent  development  of 
the  free  turbulent  mixing  layer.  Instead,  the  remnant  of  the  boundary  layer  behaves  as  a 
rotational  but  inviscid  external  flow,  and  a  new  turbulent  layer  begins  to  develop  at  the 
edge  of  this  rotational  flow. 

Several  attempts  have  been  made  to  apply  the  Korst  method  to  the  axisymmetric 
base  pressure  problem.  If  the  method  of  characteristics  is  used  to  compute  the  inviscid 
flow  adjacent  to  the  constant  pressure  portion  of  the  wake  flow,  then  the  solution 
cannot  be  extended  to  the  axis.  Therefore,  several  investigators  (e.g..  Ref.  16)  have  used 
the  assumption  that  the  flow  recompresses,  in  a  locally  planar  process,  on  an  imaginary 
sting  which  extends  from  the  base.  Mueller  (Ref.  17)  assumed  that  recompression  occurs 
either  on  a  sting  with  a  radius  of  one-half  of  the  base  radius,  or  on  a  sting  equivalent  in 
radius  to  the  experimental  radius  of  the  wake  after  recompression;  the  latter  radius  was 
determined  as  an  empirical  function  of  the  free-stream  Mach  number. 

The  simple  recompression  concepts  of  the  original  Korst  theory  are  poorly  suited  for 
the  computation  of  axisymmetric  near-wake  flows.  However,  there  are  enough  free 
parameters  in  the  axisymmetric  extensions  of  the  Korst  model  so  that  reasonably 
adequate  values  for  the  base  pressure  can  be  obtained.  But  these  models  really  are  used 
to  correlate  experimental  base  pressure  measurements  rather  than  to  predict  them. 
Moreover,  the  Korst-type  models  do  not  provide  any  information  about  the  detailed  flow 
structure  in  the  recompression  region. 

Other  methods  -  McDonald  (Ref.  18)  took  a  unique  approach  to  the  prediction  of 
the  planar  base  pressure  for  flows  which  reattach  to  a  solid  wall.  (The  centerline  of  Fig.  1 
is  replaced  by  a  wall).  McDonald  also  applied  his  analysis  to  axisymmetric  flows  in  which 
a  relatively  large  sting  extends  downstream  from  the  base  (Ref.  19).  The  base  pressure  is 
uniquely  determined  as  that  which  causes  the  reattached  flow  downstream  of 
recompression  to  become  a  fully  developed  flat  plate  boundary  layer.  Although 
McDonald  produced  reasonable  predictions  of  the  base  pressure,  his  model  can  be 
criticized  on  physical  grounds.  Numerous  experiments  have  shown  that  the  planar 
near-wake  structure  for  reattachment  to  a  plane  of  symmetry  is  not  significantly  different 
from  that  for  reattachment  to  a  wall.  Therefore,  wall  viscous  effects  cannot  be  a  major 
controlling  factor  in  the  near-wake  structure. 

Chow  (Ref.  20)  has  developed  a  two-layer  integral  treatment  of  recompression  in 
planar  flow;  in  many  aspects,  this  physically  realistic  model  combines  features  of  both 
the  Crocco-Lees  and  Korst  m.odels.  Chow  found  it  necessary  to  include  the  lateral 
pressure  gradients  in  the  recompressing  viscous  layer;  however,  this  need  to  include  lateral 
pressure  gradients  is  probably  related  to  the  assumed  geometry  of  the  dividing  streamhne. 
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This  recompression  analysis  is  applicable  only  to  the  region  upstream  of  the  RSP  (Fig.  1). 
The  predictions  of  the  Chow  model  are  discussed  in  Section  3.0. 

Recently,  Chow  and  Spring  (Ref.  21)  developed  a  method  for  predicting  the  wake 
development  downstream  of  the  RSP.  In  addition,  they  improved  the  computation  of  the 
turbulent  mixing  upstream  of  the  onset  of  recompression,  and  quite  good  predictions  of 
planar  base  pressures  were  obtained. 

Weng  (Ref.  22)  has  extended  the  Chow  recompression  model  to  axisymmetric  wake 
flows  in  which  a  sting  extends  from  the  base;  it  seems  that  this  model  should,  in 
principle,  also  apply  to  flows  without  a  sting.  In  common  with  the  models  of  Chow  and 
of  Chow  and  Spring,  the  finite  recirculation  velocities  are  neglected  in  the  flow  upstream 
of  the  recompression  region.  The  external  flow  is  computed  by  the  method  of 
characteristics  downstream  to  the  onset  of  recompression;  Prandtl-Meyer  coupling  is  used 
in  the  recompression  region.  The  initial  boundary-layer  profile  after  the  base  expansion 
process  is  taken  as  the  initial  condition  for  the  development  of  the  free  shear  layer.  Weng 
used  Prandtl  mixing  length  theory  to  compute  the  turbulent  shear  stresses  in  the  flow 
upstream  of  the  onset  of  recompression;  in  the  recompression  region,  he  used  an  arbitrary 
spatial  variation  of  the  eddy  viscosity. 

Weng  obtained  reasonably  good  predictions  of  the  effect  of  sting-base  radius  ratio  on 
the  base  pressure  for  free-stream  Mach  numbers  from  1.5  to  2.5. 

Summary  of  Earlier  Analytical  Work  -  Of  the  available  analytical  models  for  planar 
wake  flows,  those  of  Alber  (Refs.  6  and  7)  and  Chow  and  Spring  (Ref  21)  are  the  most 
perceptive.  Both  models  provide  a  prediction  of  the  tlow-field  structure  as  well  as  of  the 
base  pressure.  The  Chow  and  Spring  model  provides  the  best  prediction  of  the  effect  of 
initial  boundary  layer  on  the  base  pressure  over  the  Mach  number  range  investigated,  but 
Alber's  model  is  attractive  because  finite  reverse  flow  velocities  are  included  in  the  region 
upstream  of  the  onset  of  recompression.  Neither  of  these  models  has  been  applied  to 
flows  with  base  bleed;  indeed,  only  the  Korst-type  models  have  been  used  to  predict  the 
effects  of  base  bleed. 

Of  the  available  models  for  axisymmetric  wake  flows,  only  those  of  Strahle,  et  al. 
(Refs.  8  and  9)  and  Weng  (Ref  22)  are  based  on  a  reasonably  adequate  description  of 
the  detailed  flow  processes.  Strahle,  et  al.,  have  not  yet  included  initial  boundary-layer 
effects  in  their  model.  As  discussed  in  Section  3.0,  Weng's  method  of  computing  the 
initial  boundary-layer  effect  will  probably  be  inadequate  for  free-stream  Mach  numbers  of 
three  or  greater. 
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It  is  clear  that  even  the  recently  developed  models  for  axisymmetric  flow  do  not 
satisfy  the  criteria  established  in  Section  1.1  for  a  generalized  near-wake  analysis.  The 
void  between  the  requirements  and  the  capabilities  of  available  models  was  even  more 
pronounced  at  the  time  that  the  work  described  in  this  report  was  initiated.  Therefore,  it 
was  decided  to  undertake  the  development  of  an  analytical  model  which  would  be 
applicable  to  a  broad  range  of  supersonic  near-wake  flows  (Section  1.1). 

1.3  IMPORTANT  FEATURES  OF  A  REALISTIC  ANALYTICAL  MODEL 

When  one  considers  the  earlier  analytical  work,  as  well  as  the  available  experimental 
information;  it  becomes  clear  that  a  perceptive  model  of  the  near-wake  structure  should 
include  several  features.  First,  adequate  treatment  of  the  external  flow  in  the 
axisymmetric  case  requires  the  use  of  the  method  of  characteristics.  Second,  the  evidence 
on  quenching  of  the  turbulence  in  the  expanding  boundary  layer  suggests  that  the 
remnant  of  the  initial  boundary  layer  can  be  computed  with  the  rotational  method  of 
characteristics.  Third,  as  shown  by  Weinbaum  (Ref.  23),  the  viscous  forces  are  negligible 
compared  to  inertial  forces  during  the  turning  of  the  initial  boundary  layer,  and  the 
expansion  of  the  supersonic  portion  of  the  initial  boundary  layer  can  be  computed  with 
the  rotational  method  of  characteristics.  (A  streamtube  method  can  be  used  to  compute 
the  expansion  of  the  subsonic  portion  of  the  initial  boundary  layer).  Fourth,  a  new 
turbulent  shear  layer  should  be  considered  to  begin  after  the  corner  expansion  process.  It 
should  be  noted  that  the  concept  of  a  new  shear  layer  which  develops  downstream  of  the 
comer  expansion  has  been  suggested  by  Korst  (Ref.  24)  and  has  been  used  in  the  planar 
analysis  of  Chow  and  Spring  (Ref.  21). 

The  available  experimental  evidence  suggests  that  the  near-wake  region  between  the 
onset  of  recirculation  and  the  beginning  of  recompression  can  be  considered  to  be  a 
constant  pressure  region,  as  has  been  assumed  by  most  investigators.  However,  the 
experimental  evidence  does  not  suggest  that  the  reverse  flow  velocities  can  be  neglected, 
as  has  been  done  in  the  Korst-type  models;  the  finite  recirculation  velocities  must  be 
included. 

The  analytical  model  which  is  described  in  Section  2.0  incorporates  the  features 
which  have  been  discussed.  The  framework  of  the  analysis  is  conceptually  similar  to  the 
Crocco-Lees  approach.  However,  the  detailed  treatment  of  the  flow  processes  is  closely 
related  to  the  method  used  by  the  authors  for  ducted  two-stream  mixing  flows  (Refs.  25 
and  26). 
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2.0  DEVELOPMENT  OF  ANALYTICAL  MODEL 

The  geometry  and  nomenclature  for  the  near-wake  analysis  is  illustrated  in  Fig.  3. 
The  flow  downstream  of  the  base  is  divided  into  three  distinct  regimes:  (1)  the  "jet 
regime"  which  extends  from  the  base  plane  to  the  station  where  recirculation  begins,  (2) 
the  "wake  regime"  which  extends  from  the  onset  of  recirculation  downstream  to  the 
onset  of  recompression,  and  (3)  the  "fully  developed"  regime  which  begins  at  the  onset 
of  recompression  and  extends  downstream  into  the  fully  developed  far-wake  region. 


Figure  3.  Nomenclature  for  near-wake  analysis. 

2.1  BASIC  ASSUMPTIONS 

The  major  assumptions  used  in  the  analysis  are: 

1.  The  flow  is  steady  and  is  either  planar  or  axisymmetric. 

2.  All  gases  obey  the  perfect  gas  law. 

3.  The  integral  form  of  the  boundary-layer  equations  describes  the  flow 
between  the  axis  and  the  outer  edge  of  the  turbulent  mixing  layer,  rg. 

4.  The  turbulent  Prandtl  and  Schmidt  numbers  are  unity. 
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5.  The  free  shear  layer  which  forms  downstream  of  the  base  expansion  process 
is  fully  turbulent  (negligible  laminar  transport)  and  has  a  negligible  initial 
thickness. 

6.  The  flow  outside  of  the  free  turbulent  shear  layer  (r  >  rj)  is  supersonic, 
adiabatic,  and  inviscid. 

7.  The  upper  boundary  of  the  supersonic  external  flow  is  inviscid  and  is  either 
a  prescribed  geometry  or  has  a  prescribed  constant  pressure. 

8.  The  normalized  profiles  of  the  axial  component  of  velocity  in  the  turbulent 
shear  layer  (between  ri  and  rj)  are  similar  and  are  represented  by  a  cosine 
function. 

9.  The  element  species  and  axial  component  of  velocity  between  the  axis  and 
rj  are  one-dimensional. 

10.  The  centerline  static  pressure  in  the  jet  regime  is  an  isentropic  function 
of  Uc. 

11.  The  centerline  static  pressure  is  constant  in  the  wake  regime. 

12.  The  initial  boundary  layer  is  turbulent  and  is  characterized  by  its  thickness 
and  a  power  law  velocity  profile. 

13.  Pressure  feedback  effects  in  the  boundary  layer  upstream  of  the  base  are 
negligible. 

14.  The  expansion  of  the  initial  boundary  layer  at  the  base  plane  is  inviscid. 

15.  The  base  bleed  flow  is  uniform  and  subsonic  at  the  base  plane. 


16.  Chemical  reactions  in  the  viscous  wake  flow,  if  they  occur,  are  assumed 
to  be  in  equilibrium. 

17.  If  the  planar  flow  reattaches  to  a  wall  rather  than  to  a  plane  of  symmetry, 
then  viscous  effects  along  the  wall  are  negligible. 


The  use  of  the  boundary-layer  approximations  (assumption  3)  is  consistent  with 
Crocco-Lees  and  implies  that  (1)  there  are  negligible  lateral  pressure  gradients  in  the 
viscous  wake  flow,  and  (2)  streamwise  turbulent  transport  mechanisms  are  negligible.  The 
assumption  of  negligible  lateral  pressure  gradients  is  reasonably  valid  in  the  jet  and  wake 
regimes  because  the  curvature  of  the  streamlines  is  small  in  the  high-speed  portion  of  the 
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viscous  layer;  this  assumption  is  less  valid  in  the  recompression  region  of  the  fully 
developed  regime.  The  use  of  a  cosine  profile  in  the  shear  layer,  along  with  a  uniform 
reverse  flow  profile,  was  suggested  by  Green  (Ref.  27)  for  incompressible  wake  flows.  It 
should  be  noted  that,  even  though  the  common  term  "base  bleed"  is  used  in  this  report, 
the  momentum  of  the  stream  issuing  from  the  base  is  not  neglected;  therefore,  the 
analysis  is  applicable  to  "bleed"  velocities  in  the  high  subsonic  range. 

2.2  BASIC  EQUATIONS 

The  basic  boundary-layer  equations  are: 

Continuity 

(pur«)  +  ^  (pvr«)  =0  (1) 

ax  ar 

where 

a  =  0  for  planar  flow  and 
a  =  1  for  axisymmetric  flow 

Axial  Momentum 


pur 


du 


+  pvT 


du 


^  =  ^  (rr«)  -  r«  ^ 


dp 


(2) 


where 


r  is  the  turbulent  shear  stress 


Element  Species  Conservation 


pur 


a  dC 
dx 


+ 


pvr 


a  dC 
dr 


where  q  is  a  turbulent  species  flux. 


(3) 


Equations  (1)  through  (3)  are  integrated  (Ref.  25)  to  obtain  five  integral  equations: 
(1)  a  continuity  equation  for  the  flow  between  the  centerline  and  rj,  (2)  a  momentum 
equation  for  the  flow  between  the  centerline  and  tg,  (3)  a  momentum  equation  for  the 
flow  between  the  centerline  and  rm  (where  rm  =  rj  +  b/2),  (4)  a  species  conservation 
equation  for  the  flow  between  the  centerline  and  rj,  and  (5)  a  species  conservation 
equation  for  the  flow  between  the  centerline  and  rn,.  The  integral  equations  are: 
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Overall  Continuity 
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Overall  Momentum 
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Half-Radius  Species 
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Because  the  flow  is  one-dimensional  for  0  <  r  <  rj,  the  integral  equations  can  be 
rewritten  as  follows: 


Overall  Continuity 
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Overall  Momentum 
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Overall  Species 
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Equations  (9)  through  (13)  apply  to  all  three  regimes  of  Fig.  3.  Of  course,  q  goes  to 
zero  in  the  fully  developed  regime,  the  pressure  gradient  terms  are  zero  in  the  wake 
regime,  and  Cc  =  Cj  =  1  in  the  jet  regime. 


2.3  SHEAR  LAYER  PROFILES 


The  free  shear  layer  velocity  profile  is  given  by 


The  velocity  at  the  half-radius  surface,  I'm  (=  q  -t  b/2),  is  given  by 


u 


m 


uJ 


(14) 


(15) 


If  the  shear  layer  were  at  constant  pressure,  with  constant  values  of  Uj  and  Uc,  then 
the  element  species  concentration  profile  would  be  given  by 


C  =  A  +  B  u 


or 


(16) 


which  is  the  well  known  Crocco  integral  solution  obtained  from  Eqs.  (2)  and  (3).  In  the 
near-wake  problem,  however,  the  pressure  gradient  terms  are  negligible  only  in  the  wake 
regime.  In  addition,  Cc,  iq,  and  Us  vary  with  x,  even  in  a  region  of  constant  pressure. 
Therefore,  Eq.  (16)  cannot  be  used.  The  approach  which  has  been  taken  in  this  study  is 
to  assume  that  the  deviations  from  Eq.  (16)  can  be  expressed  as 


c 

"c 


S 


(17) 


where  K  can  vary  with  x.  It  should  be  noted  that  a  power  law  relation  similar  to  Eq. 
(17)  has  been  used  to  compensate  for  the  deviations  from  the  Crocco  relation  which  are 
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caused  by  nonunity  transport  coefficient  ratios  in  free  turbulent  flows  (Ref.  28).  Since  Cg 
is  defined  to  be  zero,  Eq.  (17)  becomes 


(18) 


2.4  COMPUTATION  OF  DENSITY  IN  THE  VISCOUS  FLOW 


For  unity  turbulent  Prandtl  and  Schmidt  numbers,  the  stagnation  enthalpy  and  the 
gas  properties  are  uniquely  related  to  C,  and  the  density  field  is  calculated  with  the 
techniques  described  in  Ref.  25.  The  stagnation  temperature,  specific  heat,  and  gas 
constant  are  specified  as  a  function  of  C,  which  varies  from  zero  (pure  outer-stream  gas) 
to  one  (pure  base  bleed  gas).  For  chemically  frozen  flow,  only  the  To,  Cp  and  R  values  at 
C  =  O  and  at  C  =  1  need  be  specified.  For  flow  with  equilibrium  chemical  reactions.  To, 
Cp  and  R  are  specified  as  functions  of  C. 

At  any  point  in  the  shear  layer,  u  is  given  by  Eq.  (14),  C  is  given  by  Eq.  (18),  and 
To,  Cp,  and  R  are  then  determined.  The  density  can  be  obtained  from  the  energy 
equation  and  the  perfect  gas  law.  Therefore,  the  density  can  be  expressed  functionally  as 

P  =  pin,  C,  p)  (19) 


2.5  TURBULENT  TRANSPORT  TERMS 

The  turbulent  shear  stress,  Tm ,  which  appears  in  Eq.  (11)  is  given  by 

dui 

r  =  p  £  — 
m  ^ m  I m 

where  e  is  the  turbulent  eddy  viscosity.  The  eddy  viscosity  is  given  by 


(20) 


£ 


(21) 


where  Rx  is  a  turbulent  Reynolds  number  based  on  the  eddy  viscosity,  the  shear  layer 
width  and  the  velocity  difference  across  the  layer;  the  model  used  to  obtain  Rj  is 
described  in  Section  2.12. 


The  equation  for  Tm  becomes,  for  the  cosine  velocity  profile. 


=  2^  Pm  hs  - 


(22) 
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The  turbulent  species  flux,  Pm ,  is  given  by 


Pm^ 


dc 


m 


or 


u^|C^K(0.5)’^"^ 


(23) 


2.6  TRANSFORMATION  OF  THE  INTEGRAL  EQUATIONS 


Jet  regime  -  Assuming  that  the  flow  conditions  along  rj  can  be  specified,  sufficient 
information  is  available  to  relate  the  terms  in  Eqs.  (9)  through  (13)  to  four  flow-field 
variables.  Since  there  are  only  four  unknowns  in  the  jet  regime,  Eq.  (13)  is  not  used.  The 
four  dependent  variables  were  selected  to  be  pc,  q,  b  and  K.  Equations  (9)  through  (12) 
can  be  transformed  into  the  following  system  (Refs.  25  and  29) 


N 


t  dx 


+  N 


dr. 

1 

2  dx 


+  N 


db 

^  dx 


+  N, 


dx 


=  N. 


(24) 


where  N  =  F  for  the  continuity  equation,  N  =  G  for  the  overall  momentum  equation,  N 
=  H  for  the  half-radius  momentum  equation,  and  N  =  I  for  the  overall  species 
conservation  equation.  The  equations  for  the  coefficients  are  presented  in  Appendix  A. 
The  numerous  auxiliary  equations  necessary  for  computation  of  the  coefficients  are 
developed  in  the  same  manner  as  described  in  Ref.  25. 

Equations  (24)  are  solved  for  the  derivatives  (dpc/dx,  drj/dx,  db/dx,  dK/dx)  by  use 
of  a  matrix  factorization  technique.  The  resulting  equations  for  the  derivatives  are 
numerically  integrated  with  a  modified  Euler  technique  (variable  step  size). 

Wake  regime  -  The  five  dependent  variables  selected  in  the  wake  regime  are  Uc,  q,  b, 
K  and  Cc.  Equations  (9)  tlmough  (13)  are  transformed  into  the  following  system: 


dr.  .1  jir  dC 

Nn—  +  N,^  +  No^  +  +  N. 


dx 


‘2  T, 


'3  d. 


'4  d. 


'5  dx 


=  N, 


(25) 


where  N  =  J  for  the  half-radius  species  equation.  The  equations  for  the  coefficients  are 
presented  in  Appendix  A. 
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Fully  developed  regime  -  The  five  dependent  variables  selected  in  the  fully  developed 
regime  are  pc,  Uc,  b,  K  and  Cc.  Equations  (9)  through  (13)  are  transformed  into  the 
following  system: 


Ni 


dx 


+  Nr 


dx 


TVT  db  dK 

+  Ng  ^  +  N4  —  +  Nc 


dx 


N. 


(26) 


The  equations  for  the  coefficients  are  presented  in  Appendix  A. 

2.7  SOLUTION  OF  THE  OUTER  FLOW  AND  COUPLING 

WITH  THE  VISCOUS  LAYER 

The  supersonic  flow  in  the  region  r  >  rg  is  computed  with  the  rotational  method  of 
characteristics.  The  "reference  plane"  technique,  in  which  the  supersonic  flow  parameters 
are  established  along  planes  normal  to  the  axis,  is  used.  The  axial  location  of  the 
reference  planes  coincides  with  the  axial  stations  used  in  the  numerical  integration  of  the 
viscous  layer  equations. 

Many  terms  which  appear  in  the  coefficients,  N,  of  the  viscous  layer  equations 
depend  on  the  flow  conditions  (and  their  derivatives  with  respect  to  x)  along  rg,  the 
interface  between  the  viscous  region  and  the  outer  supersonic  flow.  Therefore,  the 
solution  of  the  two  regions  must  be  coupled.  The  technique  used  for  coupling  the  two 
regions  is  illustrated  in  Fig.  4.  At  some  station,  x,  the  flow  in  both  regions  is  fully 
specified.  In  the  integration  technique  for  the  viscous  layer  equations,  the  step  size.  Ax,  is 
given,  as  are  the  tentative  values  of  p  and  rg  at  (x  +  Ax).  A  left  running  characteristic, 
along  which  the  flow  is  fully  defined,  is  constructed  from  the  the  last  boundary  point 
(Fig.  4).  The  right  running  characteristic  and  the  streamline,  each  of  which  passes  through 
the  new  boundary  point  and  intersects  the  left  running  characteristic,  are  constructed;  the 
equations  for  the  streamline  and  the  right  running  characteristic  are  solved  in  the  usual 
iterative  manner  to  define  the  entropy  and  flow  angle  at  the  new  boundary  point.  Then 
the  density  and  velocity  at  the  new  boundary  point  are  computed,  and  the  derivatives 
with  respect  to  x  of  these  flow  properties  along  rg  are  computed  with  an  upstream 
differencing  technique.  (These  derivatives  are  assumed  to  be  zero  at  x  =  0).  Therefore, 
sufficient  information  is  available  about  the  flow  properties  along  rg  so  that  the 
coefficients  of  the  viscous  layer  equations  at  (x  +  Ax)  can  be  evaluated. 

Within  the  framework  of  the  Euler  integration  technique,  this  boundary  point 
procedure  is  repeated,  with  different  values  of  Pc  and  rg  at  (x  +  Ax),  until  the  solution  of 
the  viscous  equations  converges.  The  solution  then  proceeds  downstream  to  a  new  axial 
station.  In  this  downstream  marching  procedure,  development  of  the  solutions  for  the 
two  flow  regions  is  fully  coupled. 
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Figure  4.  Coupling  of  the  shear  layer  and  outer-flow  solutions. 


Even  though  the  outer  flow  is  treated  with  the  rotational  method  of  characteristics, 
no  shock  waves  are  allowed  to  develop,  ffowever,  two  shock  waves  exist  even  in  simple 
near-wake  fiow  fields,  the  lip  shock  and  the  recompression  shock.  The  lip  shock,  which 
has  been  experimentally  investigated  by  Hama  (Ref.  30),  results  from  the  overturning  of 
the  low  Mach  number  portion  of  the  initial  boundary  layer  in  the  base  expansion  process 
(Ref.  23).  In  neglecting  the  lip  shock,  we  have  treated  it  as  an  isentropic  compression 
which  is  justified  since  experiments  have  shown  that  the  lip  shock  is  rather  weak,  at  least 
for  free-stream  Mach  numbers  up  to  2.0  or  3.0. 

The  series  of  compression  waves  wliich  result  from  the  recompression  process 
coalesce  to  form  the  recompression  shock  wave  (Fig.  1).  However,  the  shock  is  formed 
too  far  downstream  to  affect  the  near-wake  closure  mechanism.  Therefore,  neglecting  the 
recompression  shock  wave  influences  only  the  downstream  portion  of  the  recompression 
process. 

Inclusion  of  the  lip  and  recompression  shock  waves  in  the  present  analysis  is  quite 
feasible,  but  would  require  a  significant  increase  in  the  computation  time. 

2.8  EXPANSION  OF  THE  INITIAL  BOUNDARY  LAYER 

The  expansion  of  the  supersonic  portion  of  the  initial  boundary  layer  (Fig.  5)  is 
computed  with  the  rotational  method  of  characteristics.  Since  pressure  feedback  effects 
are  neglected,  the  supersonic  expansion  process  is  initiated  at  x  =  0  with  a  Prandtl-Meyer 
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turn  (from  pi  to  Pb)  at  the  sonic  point  in  the  boundary  layer.  Expansion  of  the  subsonic 
portion  of  the  initial  boundary  layer  is  computed  with  an  isentropic  streamtube  model. 
The  subsonic  layer  is  divided  into  ten  streamtubes,  each  of  which  expands  isentropically 
from  Pi  to  Pb.  The  flow  direction  of  the  streamtubes  after  expansion  is  taken  to  be  the 
same  as  that  of  the  initially  sonic  streamline  after  it  passes  through  the  Prandtl-Meyer 
turn.  The  flow  in  the  wall  streamtube,  after  expansion,  is  the  initial  edge  condition  (at  x 
=  0)  for  the  new  turbulent  shear  layer  which  develops  downstream  of  the  base. 


For  values  of  Pb/pi  greater  than  the  sonic  static-total-pressure  ratio,  one  or  more  of 
the  streamtubes  will  be  subsonic  after  the  expansion  process.  In  such  cases,  the  Mach 
number  of  the  subsonic  streamtubes  is  arbitrarily  adjusted  to  be  slightly  greater  than 
unity.  The  effect  of  this  adjustment  process  is  very  small  for  typical  initial 
boundary-layer  thicknesses  because  the  mass  flow  in  the  adjusted  streamtubes  is  small  and 
is  quickly  entrained  by  the  turbulent  shear  layer. 
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2.9  CHARACTERIZATION  OF  THE  INITIAL  AND  BOUNDARY  CONDITIONS 

The  turbulent  near-wake  flow  field  is  fully  determined  by  the  following  conditions, 
which  must  be  specified: 

1 .  The  upper  boundary  of  the  supersonic  external  flow, 

2.  The  base  bleed  flow  rate, 

3.  The  total  temperatures  and  gas  properties  of  the  external  flow  and  the 
bleed  flow, 

4.  Mixture  gas  properties  and  total  temperature  as  a  function  of  C,  if  chemical 
reactions  occur  (Section  2.4), 

5.  Thickness,  5,  and  velocity  profile  shape  of  the  adiabatic  initial  boundary 
layer,  and 

6.  The  profiles  at  the  base  plane  of  the  static  pressure,  total  pressure  and 
flow  angle,  between  the  edge  of  the  boundary  layer  and  the  upper  boundary. 

For  computation  of  near-wake  flows  embedded  in  external  flows  of  large  lateral 
extent  (effectively  infinite),  the  upper  boundary  is  specified  as  a  wall  parallel  to  the 
centerline  and  located  far  enough  from  the  near  wake  so  that  reflected  waves  from  the 
corner  expansion  process  do  not  interact  with  the  recompression  region. 

The  initial  boundary-layer  velocity  profile  is  assumed  to  be  given  by  a  power  law 

u/Uai  =  (y/S)" 

where  Uai  is  the  velocity  at  the  edge  of  the  boundary  layer,  6  is  the  thickness  of  the 
layer,  y  is  the  distance  from  the  wall  (Fig.  3),  and  n  is  specified. 

The  initial  boundary-layer  characteristics  for  near-wake  experiments  are  usually 
specified  in  terms  of  the  momentum  thickness  and  the  momentum  thickness  Reynolds 
number,  Re^ .  Therefore,  a  method  was  developed  to  estimate  n  as  a  function  of  M^i  and 
Reg.  Maise  and  McDonald  (Ref.  31)  showed  that  supersonic  boundary -layer  profiles  on 
an  adiabatic  flat  plate  can  be  well  represented  by  the  Coles 
law-of-the-wall/law-of-the-wake  profiles,  along  with  the  Van  Driest  transformed  velocities 
and  the  Spalding-Chi  skin  friction  values.  The  0/5  values  shown  by  Maise  and  McDonald 
as  a  function  of  Mai  and  were  matched  with  9/5  values  for  power  law  profiles  (Ref. 
32)  to  determine  n  as  a  function  of  Mgi  and  Reg  (Fig.  6).  The  relation  of  n,  Mai,  and 
0/5  is  shown  in  Fig.  7. 


26 


AEDC-TR-76-127 


The  power  law  profile,  as  well  as  the  profile  used  by  Maise  and  McDonald,  does  not 
smoothly  join  the  external  stream  velocity  profile.  The  resulting  "comer"  on  the 
total-pressure  profile  was  found  to  cause  computational  problems  in  the  viscous  layer 
solution  when  the  corner  is  entrained  into  the  shear  layer.  To  alleviate  this  problem,  the 
outer  edge  of  the  initial  boundary-layer  velocity  profile  is  smoothed  with  a  cubic 
function  which  extends  from  (rb  +  0.85)  to  (rb  +  1.25).  This  smoothing  has  a  negligible 
effect  on  the  momentum  thickness. 


Figure  7.  Boundary-layer  momentum  thickness  ratio  for 
adiabatic  flow  along  a  flat  plate. 
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2.10  TWO-LAYER  TREATMENT  OF  THE  OUTER  FLOW 

It  was  discovered  early  in  this  investigation  that  the  entropy  profile  in  the  initial 
boundary  layer  tends  to  spread  into  the  outer  flow,  which  is  usually  irrotational.  Of 
course,  the  diffusion  of  the  entropy  field  across  streamlines  is  incorrect  since  the  flow  is 
assumed  to  be  inviscid.  Apparently,  this  false  diffusion  in  the  method  of  characteristics 
flow  field  is  related  to  the  rapid  expansion  of  the  higlily  rotational  boundary-layer  flow. 
The  false  diffusion  effect  is  not  unique  to  the  present  method  of  characteristics 
formulation;  other  widely  used  method  of  characteristics  computer  programs  were  found 
to  yield  similar  results. 

To  eliminate  the  spurious  entropy  diffusion,  a  two-layer  technique  is  used  (Fig.  8). 
The  streamline  which  originates  at  the  outer  edge  of  the  inital  boundary  layer  is 
established,  and  the  total  pressure  on  this  streamline  is  specified  to  be  poa-  Therefore,  the 
entropy  field  in  the  supersonic  flow  is  divided  into  two  layers,  and  the  entropy  gradients 
from  the  initial  boundary  layer  are  confined  to  the  lower  layer. 


Left  Running  Characteristics  Which 
Originate  at  Prandtl -Meyer  Expansion 


2.11  METHOD  OF  SOLUTION 

As  in  the  Crocco-Lees  analysis,  a  saddle-point  singularity  (the  critical  point),  which 
is  intrinsic  to  the  solution,  occurs  in  the  reattachment  region.  The  behavior  of  the 
solution  is  illustrated  by  the  axial  distribution  of  centerline  velocity  (Fig.  9).  The  correct 
solution  is  that  which  passes  smoothly  through  the  singularity.  If  the  assumed  base 
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pressure  is  too  low  the  solution  will  be  supercritical  and  "blow  up",  i.e.,  in  terms  of 
Cramer's  rule,  the  determinant  of  the  coefficients  for  Eqs.  (26)  will  go  to  zero  before  the 
numerator  determinants  do.  If  the  assumed  base  pressure  is  too  high,  the  solution  will  be 
subcritical  and  is  characterized  by  the  centerline  velocity  reaching  a  maximum,  then 
decreasing  with  distance  (Fig.  9).  By  iterating  on  the  initial  conditions  (base  pressure), 
one  can  converge  on  the  set  of  initial  conditions  which  will  cause  the  flow-field  solution 
to  pass  smoothly  through  the  singularity. 


The  subcritical  solutions  are  identified  by  dUc/dx  going  to  zero.  The  supercritical 
solutions  are  identified  by  an  inflection  in  the  centerline  velocity  distribution,  i.e.,  by 
d^Uc/dx^  going  to  zero,  which  implies  that  the  determinant  is  approaching  zero.  A 
halving  technique,  with  the  bounds  established  by  the  last  subcritical  and  the  last 
supercritical  values  of  Pb,  is  used  to  converge  to  the  correct  value  of  pb.  With  a 
reasonable  first  guess  on  the  bounds  for  pb,  say  within  20  percent  of  the  correct  value, 
the  converged  value  of  Pb  can  be  identified  to  within  0.5  percent  with  approximately  ten 
iterations  of  the  flow-field  solution. 

Although  it  is  not  difficult  to  define  the  initial  conditions  with  practical  engineering 
accuracy,  establishing  the  solution  through  and  beyond  the  critical  region  is  very  difficult. 
The  nature  of  a  saddle-point  singularity  is  such  that  very  small  changes  in  the  initial 
conditions  result  in  large  changes  in  the  solution  as  the  singularity  is  approached.  Two 
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techniques  have  been  used  to  obtain  solutions  which  pass  through  the  critical  region.  The 
first  is  a  "brute  force"  approach;  the  correct  initial  condition  is  established  with  a  very 
small  tolerance,  typically  to  within  one  part  in  10^.  At  a  station  just  upstream  of  the 
critical  region,  the  centerline  velocity  is  perturbed  upward  by  a  small  amount  (typically 
one  percent)  and  the  solution  will  usually  proceed  smoothly  through  the  critical  region. 
This  procedure  has  been  found  to  work  fairly  well  for  planar  flows,  and  for  axisymmetric 
flows  with  high  external  stream  Mach  numbers.  However,  for  axisymmetric  flows  with 
moderate  external  stream  Mach  numbers  (2.0  or  3.0),  this  procedure  has  not  been 
successful,  even  when  the  initial  conditions  are  converged  to  within  one  part  in  1 0® . 

An  alternate  approach  for  developing  the  solution  through  the  critical  region  is  to 
converge  the  initial  conditions  fairly  tightly,  perhaps  to  within  one  part  in  10^  or  10^, 
then  to  extrapolate  the  centerline  velocity  distribution.  Results  at  two  axial  stations  near 
the  RSP,  which  is  well  upstream  of  the  critical  point,  are  used  to  establish  an  exponential 
extropolation  equation  for  Uc(x).  One  equation,  the  continuity  equation,  is  removed  from 
Eqs.  (26),  and  the  singular  behavior  of  the  solution  is  removed.  Experience  has  shown 
that  the  exponential  extropolation  can  be  extended  downstream,  with  good  results,  to 
the  peak  of  the  centerline  pressure  distribution  in  axisymmetric  flow. 

Results  obtained  with  both  of  the  techniques  for  developing  the  solution  through 
the  critical  region  are  presented  in  Section  3.0. 

Fortunately,  the  engineer  is  usually  interested  only  in  a  prediction  of  the  initial 
conditions  (the  base  pressure).  The  initial  conditions,  and  the  flow  field  downstream  to 
the  vicinity  of  the  RSP,  can  be  determined  with  reasonable  accuracy  without  developing 
the  solution  through  the  critical  region. 

The  occurrence  of  the  saddle-point  singularity  in  the  solution  is  a  result  of  the 
elliptic  nature  of  the  near-wake  flow  field.  That  is,  an  essentially  elliptic  problem  is 
treated  as  an  iterative  initial  value  problem,  and  the  critical  point  is  the  mechanism 
whereby  the  downstream  wake  closure  conditions  are  manifested. 

The  near-wake  analysis  has  been  programmed  for  numerical  solution  with  an  IBM 
370-165  computer.  A  typical  computation,  for  one  assumed  set  of  initial  conditions, 
which  proceeds  from  the  base  plane  to  downstream  of  the  RSP  requires  a  computation 
time  of  from  30  to  60  seconds.  The  determination  of  the  base  pressure  to  practical 
engineering  accuracy  (less  than  one  percent  uncertainty)  involves  a  computation  time  of 
five  to  ten.  minutes. 

For  near-wake  flows  in  which  there  is  no  base  bleed,  or  in  which  the  base  bleed  gas 
has  the  same  total  temperature  and  gas  properties  as  the  Outer  flow,  a  simplified  version 
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of  the  computer  program  is  used.  The  concentration  equations  are  not  solved,  and  K  and 
Cc  are  set  at  unity.  This  simplified  analysis  is  referred  to  as  the  "single  gas"  version. 

2.12  MODEL  FOR  THE  TURBULENT  REYNOLDS  NUMBER 

As  with  any  turbulent  flow  problem,  the  near-wake  analysis  requires  an  adequate 
model  for  the  turbulent  transport  properties.  Early  in  this  investigation,  an  attempt  was 
made  to  incorporate  an  integral  formulation  of  the  turbulent  kinetic  energy  (TKE) 
equation  into  the  near-wake  analysis.  This  TKE  formulation  has  been  shown  to  yield  very 
good  results  for  a  variety  of  constant  pressure-free  turbulent  flows  (Ref.  33).  But  the 
results  for  supersonic  near-wake  flows  were  very  poor;  the  turbulent  shear  stress  levels  in 
the  reattachment  region  were  grossly  underpredicted.  The  failure  of  the  TKE  formulation 
is  attributed  to  the  "extra  strain"  effects  described  by  Bradshaw  (Ref.  34).  He  points  out 
that  additional  strains  in  the  flow,  in  addition  to  the  basic  3u/3r  strain,  have  a 
disproportionately  large  effect  on  the  turbulence  structure.  In  the  recompression  region, 
the  strains  which  result  from  both  streamline  curvature  and  dilatation  (bulk  compression) 
effects  are  significant. 

Because  of  the  failure  of  the  TKE  approach,  we  decided  to  develop  an  eddy 
viscosity  model  for  the  near-wake  analysis,  even  though  the  history  of  eddy  viscosity 
models  (Ref.  35)  indicates  that  such  a  model  would  probably  have  a  limited  range  of 
applicability. 

The  eddy  viscosity  model  described  in  Section  2.5  is  the  Prandtl  model  with  the 
coefficient  expressed  as  the  inverse  of  the  turbulent  Reynolds  number,  Rj.  In  the  jet  and 
wake  regimes,  Rj  is  assumed  to  be  a  function  only  of  the  local  Mach  number,  Ms,  at  the 
supersonic  edge  of  the  shear  layer.  The  predictions  of  the  integral  TKE  method  (Ref.  33) 
agree  well  with  the  available  experiments  on  the  growth  of  planar  shear  layers  with  zero 
secondary  velocity;  the  corresponding  Rj  values  from  the  TKE  analysis  (Fig.  10)  are  well 
described  by  the  function 


R.J,  =  390  -  333  e.xp(-0.495  M  J  (27) 

which  is  used  in  the  jet  and  wake  regimes. 

The  Rt  model  used  in  the  fully  developed  regime  is  based  on  the  experimental 
results  of  Mueller  (Ref.  36)  and  Tani  (Ref  37)  for  incompressible  reattaching  planar 
flows.  Their  experimental  configurations  are  shown  in  Fig.  1 1 .  The  experimental  velocity 
profiles  were  approximately  fitted  with  a  cosine  function  to  estimate  the  width  of  the 
shear  layer  and  the  effective  wall  slip  velocity.  The  results,  along  with  the  measured  shear 
stress  at  the  center  of  the  velocity  profile,  were  used  to  compute  Rj.  Two  characteristic 
Rx  values  are  of  interest,  the  value  at  the  beginning  of  reattachment,  Rjp  and  the  value 
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at  the  RSP,  Rxi.  (The  beginning  of  reattachment  is  defined  as  the  axial  station  where  the 
wall  static  pressure  is  equal  to  the  base  pressure.)  The  values  for  Rjp  and  Rxr  are  shown 
below. 


Experiment  Rxp  Rxr 


Mueller,  h  = 

0.75  in. 

276 

119 

Mueller,  h  = 

0.50  in. 

216 

106 

Mueller,  h  = 

0.25  in. 

219 

101 

Tani,  h  =  2 

cm 

295 

108 

Because  of  the  consistent  values  for  Rxr,  from  100  to  120,  we  speculated  that  Rjr 
would  be  a  good  parameter  to  characterize  the  level  of  the  eddy  viscosity  in  the 


reattachment  region. 


Figure  10.  Turbulent  Reynolds  number  for  planar  shear  layers  with  no  secondary  flow. 
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Figure  11.  Configurations  for  low-speed  separating  and  reattaching 
flow  experiments. 

From  the  experimental  results  of  Mueller  and  of  Tani,  it  was  found  that  the 
distribution  of  Rx  in  the  reattachment  region  could  be  taken  to  be  a  unique  function  of 
Uc/us;  Us  is  the  local  velocity  at  the  edge  of  the  viscous  layer  and  Uc  is  the  effective  wall 
slip  velocity.  Results  are  shown  in  Fig.  12  for  the  region  between  the  beginning  of 
reattachment  and  the  RSP.  The  data  are  well  represented  by 

(R.J.  -  Rxr^^^^Tp  “  +  0.4  11"^  (28) 

where  u  =  (Uc/us)/(uc/us)p.  Downstream  of  the  RSP,  the  data  (Fig.  13)  are  well  represented 
by 

Rx/Rpr  =  exp  (-2.5  u^/uj  (29) 

In  the  near-wake  analysis,  Eqs.  (28)  and  (29)  are  used  to  express  the  axial  variation 
of  Rx,  and  Rxp  is  given  by  Eq.  (27).  Therefore,  the  level  of  the  eddy  viscosity 
throughout  the  fully  developed  regime  is  scaled  by  Rxr- 
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Figure  13.  Turbulent  Reynolds  number  downstream  of  the  RSP. 


A  series  of  computational  experiments  was  carried  out  to  establish  the  values  of 
Rxr  required  to  give  a  reasonably  accurate  prediction  of  experimental  base  pressures  in 
supersonic  flow  (say  to  within  five  percent).  The  following  matrix  of  test  cases  was 
selected: 

planar,  M  ,  =  2.0,  6/h  =  0.016 
ai 

planar,  M  =  2  0>  d/h  =  0.062 
a  1 

planar,  jV]^j  =  3.0, d/h  =  0.013 

Axisymmetric,  =  2.0,  0/h  =  0.024 
Axisymmetric,  M^j  ~  2.0,  0/h  =  0.053 
Axisymmetric,  =  4.0,  9/h  =  0.011 

Typical  effects  of  Rxr  on  the  base  pressure,  for  planar  flow  at  Mai  “  2,  are  shown  in 
Fig.  14.  As  would  be  expected,  an  increase  in  Rxr  (a  decrease  in  the  turbulent  shear 
stress  level)  causes  an  increase  in  the  base  pressure. 
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Turbulent  Reynolds  Number  at  RSP,  Rj^ 

Figure  14.  Influence  of  the  turbulent  Reynolds  number  at  the  RSP  on  the 
predicted  base  pressure. 


From  the  matrix  of  base  pressure  test  cases,  it  was  found  that  Rxt,  must  decrease  as 
the  initial  boundary-layer  thickness  is  increased  and  as  Mai  is  increased.  In  addition,  it 
was  found  that  accurate  prediction  of  axisymmetric  flows  requires  a  larger  value  of  Rxr 
than  is  required  for  planar  flows.  Since  Bradshaw's  curvature  and  dilatation  strains  are 
both  related  to  a  pressure  gradient  parameter,  we  decided  to  correlate  Rxr  with  a 
parameter  proportional  to  the  pressure  gradient  at  the  beginning  of  re  compression.  After 
some  trial  and  error,  the  pressure  gradient  parameter  was  selected  to  be 

where  5p  is  the  distance  between  the  axis  and  the  streamline  which  originates  at  the  edge, 
of  the  initial  boundary  layer.  The  results  of  the  correlation  of  Rxr  with  P  are  shown  in 
Fig.  15.  Both  curves  wliich  have  been  fitted  to  the  required  Rxr  values  can  be  expressed 
by  the  following  single  equation: 

Rnp^  =  60  +  15a  +  (40  +  135a)  exp  (-0.024  P^'^)  (31) 
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Figure  15.  Turbulent  Reynolds  number  at  the  RSP  as  a  function  of  the  pressure 
gradient  parameter. 


It  is  interesting  that  the  required  Rti  values  for  planar  flows  (Fig.  15)  at  low  values 
of  P  are  consistent  with  the  low-speed  experimental  results  of  Mueller  and  of  Tani. 

It  will  be  shown  in  Section  3.0  that  the  Rj  model  developed  in  this  section  yields 
good  results  for  planar  near-wake  flows;  in  addition  to  the  base  pressure,  the  flow 
structure  is  predicted  with  adequate  accuracy.  However,  the  Rj  model  has  been  found  to 
yield  only  fair  results  for  axisymmetric  flows  with  M^i  less  than  about  1.7;  the 
predictions  are  fairly  adequate  at  higher  Mach  numbers. 


3.0  DISCUSSION  OF  RESULTS 


In  this  section,  the  predictions  of  the  near-wake  model  are  compared  with 
experiments!  results  for  plsnsr  snd  sxisymmetric  flows  of  tlie  type  shown  in  Fig.  i.e., 
the  supersonic  external  flow  upstream  of  the  base  and  outside  of  the  wall  boundary  layer 
is  uniform  and  parallel  to  the  axis  and  is  effectively  infinite  in  lateral  extent. 
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3.1  CONCENTRATION  FIELD  SOLUTION  PROBLEMS 

When  attempts  were  made  to  compute  near-wake  flows  for  which  the  concentration 
equations  must  be  solved,  a  singularity  in  the  wake  regime  solution  was  encountered. 
Downstream  of  the  onset  of  recirculation,  the  predicted  centerline  concentration,  Co, 
decays  with  axial  distance.  At  some  axial  station  upstream  of  the  onset  of  recompression, 
the  system  of  equations  becomes  singular;  i.e.,  the  determinant  of  the  system  of 
equations  [Eq.  (25)]  goes  to  zero  and  the  solution  "blows  up".  This  singularity  is  not 
related  to  the  Crocco-Lees  saddle-point  singularity  which  occurs  in  the  recompression 
region.  Indeed,  the  new  singularity  does  not  exhibit  saddle-point  behavior;  the  solutions 
for  all  of  the  derivatives  (dCc/dx,  etc.)  remain  smooth  and  monotonic  as  the  singularity  is 
approached. 

The  wake  regime  singularity  does  not  appear  to  have  physical  significance,  but  is 
probably  related  to  the  particular  formulation  of  the  half-radius  species  equation  or  of 
the  concentration  profile  shape  function.  The  singularity  also  occurs  when  the 
concentration  equations  are  solved  for  a  single  gas  base  bleed  problem,  even  though  no 
wake  regime  singularity  occurs  when  the  single  gas  version  of  the  analysis  is  used  for  the 
same  problem. 

Strahle  and  Mehta  (Ref.  8)  encountered  an  extraneous  singularity  in  their  analysis  of 
near-wake  flows  without  base  bleed.  They  were  able  to  remove  the  singular  behavior  by 
changing  their  velocity  profile  functions  in  the  region  upstream  of  recompression  (Ref. 
9).  Various  modifications  of  the  concentration  profile  function  in  the  present  analysis 
have  been  tried,  but  to  no  avail;  the  wake  regime  singularity  occurs  for  both  single  gas 
flows  and  for  flows  in  which  the  base  bleed  gas  is  different  from  the  outer-stream  gas. 
An  obvious  next  step  in  attempting  to  solve  tliis  concentration  field  solution  problem  is 
to  alter  the  formulation  of  the  half-radius  species  equation. 

Because  of  the  unresolved  problems  in  solving  the  concentration  equations,  only 
results  from  the  single  gas  version  of  the  near-wake  model  are  presented  in  this  report. 

3.2  PLANAR  BASE  PRESSURE 

The  effect  of  initial  boundary-layer  momentum  thickness  on  base  pressure  in  planar 
flow  is  shown  in  Fig.  16  for  Mai  ~  1-5  and  in  Fig.  17  for  Mgi  =  2.0.  The  scatter  shown 
in  the  experimental  results  is  typical  for  planar  flows  and  is  attributed  to  sidewall  effects 
which  vary  from  experiment  to  experiment.  At  both  free-stream  Mach  numbers,  the 
present  analysis  predicts  the  experimental  base  pressures  very  well  for  0/h  less  than  0.03 
or  0.04.  However,  the  analysis  overpredicts  the  base  pressure  for  higher  values  of  0/h;  the 
predicted  values  of  Pb/pi  are  about  ten  percent  too  high  at  0/h  =  0.1.  In  agreement  with 
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experimental  observations  (Ref.  38)  at  free-stream  Mach  numbers  of  1.5  and  2,0,  the 
near-wake  model  predicts  a  negligible  influence  of  Re^  on  the  base  pressure  for  a  given 
value  of  0/h.  (Computations  were  made  for  Re^  from  10^  to  10^). 

The  trend  of  Pb/Pi  versus  0/h  predicted  by  the  present  analysis  is  similar  to  that 
predicted  by  Alber  (Figs.  16  and  17);  however,  the  present  analysis  provides  significantly 
better  predictions  of  the  level  of  base  pressure.  The  Chow  theory  (Ref.  20)  predicts  an 
incorrect  trend  of  Pb/Pi  versus  0/h;  however,  the  Chow  and  Spring  extension  of  this 
theory  provides  very  good  predictions  of  the  base  pressure  at  Mach  numbers  of  1.5  and 
2.0  (Figs.  16  and  17). 


Figure  16.  Effect  of  the  initial  boundary  layer  on  the  base  pressure  in  planar 
flow  (Mai  =  1.5). 

Experimental  base  pressure  results  for  Mai  =3.0  (Fig.  18)  show  much  more  scatter 
than  was  observed  at  the  lower  Mach  numbers.  In  agreement  with  the  trend  of  the 
experimental  results,  the  present  analysis  predicts  a  much  larger  effect  of  0/h  on  the  base 
pressure  at  Mai  =  3.0  than  at  Mai  “  2.0.  In  addition,  a  significant  effect  of  Re^  on  the 
base  pressure  is  predicted. 

The  effect  of  free-stream  Mach  number  on  the  base  pressure  is  shown  in  Fig.  19a  for 
0/h  =  0.02  and  in  Fig.  19b  for  0/h  -  0.04.  In  view  of  the  scatter  in  planar  base  pressure 
data,  the  predictions  shown  in  Fig.  19  are  very  satisfactory. 
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Most  near-wake  flows  of  engineering  interest  involve  initial  boundary  layers  of 
moderate  thickness,  with  0/h  values  less  than  0.05.  The  base  pressure  predictions  of  the 
present  analysis  for  this  range  of  0/h  are  shown  in  Fig.  20  for  Mai  from  1.2  to  4.0.  Also 
shown  in  Fig.  20  are  the  predictions  of  the  Korst  theory  for  no  initial  boundary  layer. 
Because  of  a  fortunate  combination  of  physically  unrealistic  assumptions,  this  simplest 
version  of  the  Korst  theory  provides  fairly  satisfactory  predictions  of  planar  base  pressure 
in  the  range  of  initial  boundary  layers  most  often  encountered  in  engineering 
applications. 


Figure  20.  Predicted  base  pressure  for  planar  flow. 


3.3  PLANAR  FLOW-FIELD  STRUCTURE 

Predicted  axial  distributions  of  the  total  pressure  at  the  supersonic  edge  of  the 
turbulent  shear  layer  are  shown  in  Fig.  21  for  three  typical  planar  wake  flows.  A  value  of 
Pos/Poa  less  than  unity  indicates  that  the  shear  layer  is  growing  into  the  rotational 
remnant  of  the  initial  boundary  layer,  and  the  point  at  which  pos/Poa  becomes  unity  is 
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the  point  where  the  initial  boundary-layer  flow  has  been  completely  entrained.  The 
results  shown  in  Fig.  21  indicate  that,  even  for  relatively  thin  initial  boundary  layers  (d/h 
=  0.013  -  0.014),  the  entire  viscous  near-wake  structure  downstream  to  the  critical  point 
is  immersed  in  a  rotational  external  flow. 


Figure  21.  Total  pressure  at  the  edge  of  the  shear  layer  for  planar  flow. 


Predicted  and  experimental  centerline  (wall)  static  pressure  distributions  are  shown 
in  Fig.  22  for  Mai  “  2.0.  The  distance  to  the  beginning  of  recompression  is  somewhat 
overpredicted  by  the  present  theory;  however,  the  complete  predicted  pressure 
distribution  is  considered  to  be  in  satisfactory  agreement  with  the  experiment,  as  is  the 
predicted  location  of  the  RSP.  The  critical  point  in  this  flow  occurs  at  a  centerline  Mach 
number  of  0.64. 

The  predicted  pressure  distribution  of  Chow  and  Spring,  which  is  shown  in  Fig.  22, 
is  fairly  satisfactory.  Their  analysis  is  earned  out  in  two  parts:  (1)  for  the  recompression 
process  up  to  the  RSP  and  (2)  for  the  redeveloping  flow  downstream  of  the  RSP.  As 
shown,  the  two  segments  of  the  predicted  pressure  distribution  do  not  join  smoothly  at 
the  RSP. 
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The  experimental  centerline  (wall)  static  pressure  distribution  obtained  by  Reda  and 
Page  (Ref.  40)  for  Mgi  =  3.05  is  shown  in  Fig.  23,  along  with  the  predicted  pressure 
distribution.  The  overprediction  of  the  distance  to  the  onset  of  recompression  is 
somewhat  greater  than  at  Mai  =2.0;  however,  the  axial  distance  to  the  experimental  RSP 
is  predicted  to  within  ten  percent.  The  centerline  Mach  number  distribution  for  the  Reda 
and  Page  flow  is  shown  in  Fig.  24.  As  a  result  of  neglecting  the  streamwise  transport 
terms  in  the  viscous  wake  flow,  the  predicted  centerline  Mach  number  distribution  always 
exhibits  the  sharp  minimum  at  the  beginning  of  recirculation.  However,  the  solution 
seems  to  recover  from  this  localized  unrealistic  behavior,  and  the  downstream  portion  of 
the  predicted  M^  distribution  is  reasonably  good.  Detailed  lateral  profiles  of  the  pitot 
pressure  in  the  Reda  and  Page  flow  are  shown  in  Fig.  25  for  four  axial  stations;  the  first 
two  stations  are  upstream  of  recompression,  and  the  last  two  are  downstream  of  the 
point  where  the  wake  flow  is  completely  supersonic.  The  overall  shapes  of  the  predicted 
profiles  agree  fairly  well  with  the  experiment,  but  it  is  evident  that  the  experimental 
shear  layer  growth  rate  in  the  wake  regime  is  underpredicted.  Part  of  the  deviation 
between  prediction  and  experiment  at  the  two  downstream  stations  may  be  attributed  to 
the  presence  of  the  strong  recompression  shock  wave,  which  is  treated  as  an  isentropic 
compression  in  the  analysis. 

The  "brute  force"  technique  described  in  Section  2.1 1  was  used  to  obtain  the  planar 
flow-field  solutions  which  extend  through  the  critical  region. 


Figure  22.  Centerline  pressure  distribution  for  planar  wake  flow  (Mai  =  2.0). 


44 


Mach 


Figure  24.  Centerline 


number  distribution  for 


wake  flow  (Mgi  =  3.05). 


Distance  from  Centerline,  r/h 


AEDC-TR-76-127 


3.4  AXISYMMETRfC  BASE  PRESSURE 

In  order  to  avoid  possible  support  interference  effects  on  the  base  pressure  data, 
results  from  experiments  without  base  stings  (free  flight  or  annular ,  tunnel  experiments) 
have  been  used  for  comparison  with  the  present  analytical  model.  The  effect  of  initial 
boundary-layer  momentum  thickness  on  the  base  pressure  is  shown  in  Fig.  26  for  Mai  = 
2.0.  As  was  found  for  planar  flows,  the  predicted  axisymmetric  base  pressure  at  Mai  ~ 
2.0  is  insensitive  to  changes  in  Reg  in  the  range  from  10^  to  10^.  Although  the 
experimental  base  pressures  are  predicted  with  less  than  ten-percent  error,  the  analytical 
results  exhibit  a  peculiar  oscillation  about  the  expected  results,  which  are  indicated  by 
the  dashed  line.  The  oscillation  is  attributed  to  the  Rx  model  used  in  the  reattachment 
region,  which  is  overly  sensitive  to  small  changes  in  the  flow  field.  The  Rx  model  yields 
realistic  results  for  higher  Mach  numbers,  say  3.0  or  4.0,  but  the  results  at  Mach  numbers 
less  than  2.0  exhibit  more  oscillation  than  shown  in  Fig.  26. 


Figure  26.  Effect  of  the  initial  boundary  layer  on  the  base  pressure  in 
axisymmetric  flow  (Mai  =  2.0). 
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The  effect  of  initial  boundary-layer  momentum  thickness  on  the  base  pressure  is 
shown  in  Fig.  27  for  Mai  =  3.9;  the  prediction  is  reasonably  satisfactory. 


Figure  27.  Effect  of  the  initial  boundary  layer  on  the  base  pressure  in 
axisymmetric  flow  (Mai  =  3.9). 


The  effect  of  free-stream  Mach  number  on  the  axisymmetric  base  pressure  is  shown 
in  Fig.  28.  In  view  of  the  considerable  scatter  in  the  experimental  data,  the  predictions 
for  0/rb  =  0.01  and  0/rb  =  0.05  are  considered  to  be  satisfactory  for  Mai  greater  than 
2.0.  However,  at  lower  Mach  numbers  the  peculiarities  caused  by  the  Rj  model  are 
evident.  The  prediction  for  0/rb  =  0,01  is  realistic  down  to  a  Mach  number  of  1.5,  below 
which  the  the  base  pressure  ratio  is  predicted  to  be  constant.  A  similar  flattening  of  the 
Pb/pi  curve  for  0/ri,  =  0.05  occurs  at  Mach  numbers  less  than  2.0.  Because  of  the 
unrealistic  behavior  of  the  recompression  region  turbulent  transport  model,  the  near-wake 
analysis  in  its  present  form  is  not  recommended  for  use  at  free-stream  Mach  numbers  less 
than  about  1.7. 
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Figure  28.  Mach  number  effect  on  the  base  pressure  for  axisymmetric  flow. 


3.5  EFFECT  OF  BASE  BLEED  ON  AXISYMMETRIC  BASE  PRESSURE 

Results  of  the  axisymmetric  base  bleed  experiments  of  Sirieix,  et  al.  (Ref.  43),  are 
shown  in  Fig.  29,  along  with  the  predictions  of  the  near-wake  analysis.  The  experimental 
trend  of  increasing  base  pressure  with  increasing  base  bleed  flow  is  well  predicted;  in 
addition,  the  maximum  difference  between  the  predicted  and  experimental  base  pressures 
is  less  than  eight  percent. 
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Figure  29.  Effect  of  base  bleed  on  axisymmetric  base  pressure. 


3.6  AXISYMMETRIC  FLOW-FIELD  STRUCTURE 

Predicted  axial  distributions  of  the  total  pressure  at  the  supersonic  edge  of  the 
turbulent  shear  layer  are  shown  in  Fig.  30  for  three  typical  axisymmetric  wake  flows.  As 
was  found  for  planar  flows,  most  or  all  of  the  viscous  near-wake  structure  is  imbedded  in 
a  rotational  external  flow.  The  abrupt  increase  in  the  slope  of  the  pos  versus  x  curves  at 
the  onset  of  recompression,  which  is  much  more  pronounced  for  axisymmetric  flow  than 
for  planar  flow,  indicates  a  sudden  increase  in  the  entrainment  rate  of  the  shear  layer. 
The  mass  flow  in  the  shear  layer  increases  rapidly  with  distance  in  the  recompression 
region,  even  though  the  thickness  of  the  shear  layer,  b,  initially  decreases  in  this  region; 
this  phenomenon  is  related  to  the  rapid  change  in  the  shear  layer  velocity  profile  during 
recompression. 
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Figure  30.  Total  pressure  at  the  edge  of  the  shear  layer  for  axisymmetric  flow. 

The  experimental  centerline  static  pressure  distribution  measured  by  Seiling,  et  al. 
(Ref.  47),  for  =  3.91  is  shown  in  Fig.  31,  along  with  the  predictions  of  the 
near-wake  analysis.  The  solution  through  the  critical  region  was  obtained  with  the  "brute 
force"  technique  described  in  Section  2.11.  The  distance  to  the  onset  of  recompression  is 
overpredicted  as  is  the  distance  to  the  RSP.  The  characteristic  overshoot  of  the  static 
pressure  is  predicted. 

The  length  scale  of  the  near  wake,  as  indicated  by  the  distance  to  the  onset  of 
recompression,  is  strongly  influenced  in  high  Mach  number  flows  by  the  rotational 
boundary-layer  flow  at  the  base  plane.  To  illustrate  this  point,  a  computation  of  the  flow 
field  of  Seiling,  et  al.,  was  made  with  the  assumption  that  there  was  no  intial  boundary 
layer,  i.e,  that  the  turning  of  the  flow  at  the  base  plane  corresponds  to  a  Prandtl-Meyer 
expansion  of  the  inviscid  free  stream  through  the  pressure  ratio  Pb/Pi-  The  results,  which 
are  presented  in  Fig.  31,  show  that  the  distance  to  the  onset  of  recompression  is  nearly 
twice  as  large  as  that  predicted  when  the  initial  boundary  layer  is  included. 
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Axial  Distance  from  Base,  xlr^ 

Figure  31.  Centerline  pressure  distribution  for  axisymmetric  wake  flow  (Mai  =  3.91). 

Detailed  radial  profiles  of  pitot  pressure,  which  were  measured  by  Seiling,  et  al.,  are 
shown  in  Fig.  32.  The  predicted  profiles  agree  well  with  the  experiment,  although  the 
experimental  shear  layer  growth  rate  in  the  wake  regime  is  underpredicted;  a  similar 
result  was  found  for  planar  flow  (Fig.  25).  Based  on  the  results  shown  in  Fig.  32,  as  well 
as  those  shown  in  Fig.  25  for  planar  flow,  it  does  not  appear  that  the  neglect  of  the 
finite-strength  lip  shock  is  a  serious  deficiency  in  the  present  analysis,  even  for 
free-stream  Mach  numbers  up  to  4.0. 

The  experimental  centerline  pressure  distribution  measured  by  Sirieix,  et  al.  (Ref. 
43),  for  Mai  ==  1.92  is  shown  in  Fig.  33.  For  Mach  numbers  in  the  neighborhood  of  2.0, 
the  "brute  force"  technique  for  obtaining  the  solution  tlirough  the  critical  region  is 
unsuccessful,  so  the  extrapolation  technique  described  in  Section  2.11  was  used  to  obtain 
the  solution  shown  in  Fig.  33.  The  distance  to  the  onset  of  recompression  is  somewhat 
overpredicted,  and  the  distance  to  the  RSP  is  overpredicted  by  25  percent.  An  unusual 
inflection  in  the  centerline  pressure  distribution  is  predicted  in  the  region  of  the  RSP,  but 
the  characteristic  overshoot  in  the  pressure  distribution  is  preserved  in  the  solution. 
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The  centerline  Mach  number  distribution  for  the  flow  of  Sirieix,  et  ah,  is  shown  in 
Fig.  34.  As  was  shown  for  planar  flow,  the  predicted  Mach  number  distribution  has  a 
sharp  minimum;  this  unrealistic  behavior  is  localized,  and  the  slope  of  the  predicted 
centerline  Mach  number  distribution  agrees  well  with  the  experiment  at  the  RSP, 
although  the  entire  distribution  is  shifted  downstream  from  the  experimental  distribution 
by  the  same  distance  as  the  RSP  shift. 

The  experimental  centerline  pressure  distribution  measured  by  Reid  and  Hastings  for 
Mai  =  2.03  is  shown  in  Fig.  35.  The  inflection  shown  in  Fig.  33  does  not  appear  in  this 
solution,  which  was  also  obtained  with  the  centerline  velocity  extrapolation  technique. 

Based  on  the  results  shown  in  Figs.  33  through  35,  it  can  be  concluded  that  the 
extrapolation  technique  for  obtaining  solutions  through  the  critical  region  is  reasonably 
successful. 


Figure  34.  Centerline  Mach  number  distribution  for 
axisymmetric  wake  flow  (Mgi  =  1.92). 
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Axial  Distance  from  Base,  x/rjj 

Figure  35.  Centerline  pressure  distribution  for  axisymmetric  wake  flow  (Mgi  =  2.03). 

3.7  DIVIDING  STREAMLINE  PROPERTIES 

Although  dividing  streamline  concepts  play  no  role  in  the  present  near-wake  analysis, 
the  location  of  this  streamline  and  its  velocity  are  computed.  Predicted  distributions  of 
total  pressure  on  the  dividing  streamline  are  shown  in  Fig.  36  for  three  typical  planar 
flows.  Somewhat  surprisingly,  the  predicted  dividing  streamline  total  pressure  remains 
nearly  constant  from  the  base  plane  to  the  RSP. 

The  simple  Korst  theory  for  no  initial  boundary  layer  yields  base  pressure  results 
which  are  similar  to  those  of  the  present  analysis  for  0/h  from  0.01  to  0.05  (Fig.  20). 
Korst 's  theory  is  based  on  three  fundamental  assumptions:  (1)  an  isentropic 
recompression  along  the  dividing  streamline,  (2)  stagnation  of  the  dividing  streamline  at 
the  peak  of  the  recompression  pressure  rise,  and  (3)  negligible  reverse  flow  velocities. 
According  to  the  results  shown  in  Fig.  36,  the  first  assumption  is  a  good  one  for  planar 
flows.  It  has  long  been  known  that  the  dividing  streamline  does  not  stagnate  at  the  peak 
of  the  pressure  rise,  i.e.,  the  RSP  occurs  at  a  station  where  the  pressure  is  much  less  than 
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the  peak  value  (Figs.  22  and  23).  But  Korst's  assumption  of  negligible  reverse  flow 
velocities  leads  to  overprediction  of  the  dividing  streamline  total  pressure  (Fig.  36).  This 
high  total  pressure  in  Korst's  analysis,  combined  with  the  unrealistically  high  pressure  at 
which  the  dividing  streamline  is  assumed  to  stagnate,  leads  to  base  pressure  predictions 
which  are  about  correct. 


Axial  Distance  from  Base,  x/h 

Figure  36.  Total  pressure  along  the  dividing  streamline  for  planar  flows. 


Predicted  distributions  of  the  total  pressure  on  the  dividing  streamline  are  shown  in 
Fig.  37  for  three  typical  axisymmetric  wake  flows.  In  contrast  to  the  results  for  planar 
flows,  the  dividing  streamline  total  pressure  is  predicted  to  increase  significantly  during 
the  recompression  process.  However,  the  quality  of  the  detailed  predictions  of  the  present 
model  for  axisymmetric  reattaching  flows  is  not  considered  to  be  good  enough  so  that 
strong  conclusions  can  be  drawn  about  the  validity  of  isentropic  recompression  of  the 
dividing  streamline. 

One  point  is  clear  from  the  results  discussed  in  this  section:  any  wake  analysis  which 
is  based  on  dividing  streamline  concepts  (e.g.,  Ref.  21)  should  include  the  effect  of  finite 
reverse  flow  velocities  on  the  dividing  streamline  properties. 
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4.0  CONCLUDING  REMARKS 

The  correlations  of  the  present  near-wake  analysis  with  a  variety  of  planar  and 
axisymmetric  experiments  indicate  that  the  overall  framework  of  the  analysis  is  physically 
realistic.  Because  both  the  detailed  flow-field  structure  and  the  base  pressure  can  be 
predicted  reasonably  well,  it  can  be  concluded  that  the  arbitrariness  in  the  selected 
turbulent  transport  model  is  not  used  to  mask  the  effect  of  physically  unrealistic  aspects 
of  the  analysis.  The  method  of  computing  the  inviscid  remnant  of  the  initial  boundary 
layer  is  considered  to  be  a  very  important  feature  of  the  analysis,  as  is  the  inclusion  of 
finite  reverse  flow  velocities. 

In  any  analysis  of  turbulent  flow,  the  selection  of  an  adequate  turbulent  transport 
model  is  a  major  problem.  The  transport  model  used  in  the  present  analysis  yields 
acceptable  base  pressure  results  for  planar  flows,  but  yields  inconsistent  results  for 
axisymmetric  flows  with  free-stream  Mach  numbers  less  than  about  2.0.  The  concept  of 
scahng  the  recompression  region  turbulent  transport  properties  with  the  value  of  the 
turbulent  Reynolds  number  at  the  rear  stagnation  point  seems  to  be  satisfactory.  But  the 
model  used  to  relate  this  turbulent  Reynolds  number  to  the  flow-field  properties  must  be 
improved  if  adequate  predictions  of  axisymmetric  flows  are  to  be  obtained  over  the 
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entire  supersonic  Mach  number  range  of  practical  interest.  The  correlations  of  the 
analytical  results  with  the  flow-field  experiments  of  Reda  and  Page  and  of  Selling,  et  ah, 
show  that  the  turbulent  mixing  rate  is  somewhat  underpredicted  in  the  region  upstream 
of  recompression  and  indicate  that  the  Rx  model  for  tliis  region  should  be  improved. 
However,  based  on  limited  numerical  experimentation,  it  is  not  expected  that  an 
improved  Rx  model  for  the  wake  regime  would  greatly  change  the  predicted  base 
pressures;  the  predicted  base  pressure  is  much  more  strongly  dependent  on  the  turbulent 
transport  properties  in  the  recompression  region  than  on  those  properties  in  the  wake 
regime. 

Because  of  the  unresolved  problems  in  the  solution  of  the  species  conservation 
equations,  the  near-wake  model  as  presently  formulated  cannot  be  used  for  flows  in 
which  the  base  bleed  gas  is  different  from  the  external  stream  gas.  Many  practical  wake 
flows  involve  two  gases;  therefore,  a  major  emphasis  in  future  work  on  the  near-wake 
model  should  be  placed  on  developing  a  satisfactory  technique  for  solving  the  species 
conservation  equations. 

Although  the  overall  framework  of  the  near-wake  analysis  seems  to  be  physically 
realistic,  the  behavior  of  the  solutions  in  the  vicinity  of  the  onset  of  recompression  is  not 
satisfactory.  The  length  of  the  constant  pressure  region  is  overpredicted  and  an  unrealistic 
sharp  minimum  in  the  centerline  velocity  is  predicted.  Both  of  these  characteristics  of  the 
solution  can  probably  be  avoided  by  abandoning  the  assumption  that  the  pressure  is 
constant  in  the  wake  regime.  In  order  to  solve  for  the  pressure  distribution,  however,  an 
additional  equation  is  required.  A  reasonable  approach  would  be  to  add  an  axial 
momentum  equation  for  the  flow  along  the  axis;  this  equation  would  include  an 
approximate  term  to  represent  the  streamwise  turbulent  transport  of  momentum. 

Even  though  the  present  formulation  of  the  near-wake  model  is  deficient  in  several 
respects,  the  model  provides  a  significant  improvement  in  the  capability  to  analyze 
supersonic  turbulent  near-wake  flows.  Indeed,  the  present  model  for  axisymmetric  flows 
without  base  stings  is  the  first  to  realistically  include  the  effect  of  the  initial 
boundary-layer  thickness  and  shape  and  the  effects  of  base  bleed. 
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APf>ENDIX  A 

COEFFICIENTS  OF  THE  INTEGRAL  EQUATIONS 


The  following  parameters  are  defined: 


y  =  (r  -  rj)/b 


Tg  =  rj  +  b 


^  dr  =  ab^ydy  +  r?b^  ;  a  =  0  or  1 

»  i 

(5  ,  2\ 

^2g  ^  ^  ^ 

^3g  =  (puC) 


The  jet  regime  coefficients  (F,  G,  H,  and  I),  where  g  =  p,  rj,  b,  K,  and  x,  are  (the 
subscript  c  on  Pc  has  been  omitted  for  convenience): 

.a+l 


Fi  =  ab^/  Sipydy  +  r«b /s^^dy  ,  ^  (p.ud 


Fg  =  ah^f  Sj^  ydy  +  r?b  /  S^^  dy 


Fg  =  ah^  f  Sgj^ydy  +  rfbf  Sjj^dy 


F4  =  ah^  f  Sjjj-ydy  +  r?b  /  S^^dy 


F5  =-ab2/  S^^ydy  +  r?b  /  Sj^dy  -  p^v 


r« 
s  s 


j,a+l  ^a+l 


Gi  =  ab  /  S2pydy  +  r?b  J  S2pdy  + 


'  d  .,2 


(P,uf) 


-(„+!)  -  («  +  i)  ^ 


G2  =  ah^  f  S2,.ydy  +  r?b  /  82^,  dy 


Gg  =  ab^  f  S2i3ydy  +  r?b  /  S2i,dy 
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G4  =  ab^  /  S2j.ydy  +  rfb  /  S 

2Y.^y 

O  o 

G5  =  -  ab2/^S2^ydy  -  r?b/^S2^dy  - 

o  O 


‘/a 


'/a 


Hi  =  ab  /  S2  ydy  +  83  dy  -  a  u^b  /  Sjpydy 

0*0  o 

a+1  .a+1 

„  ‘  d  /  2n  i  d  ,  , 

-  +  teTl)  -  V(^TT) 

H2  =  ah^ 82,  ydy  +  rfb/'82,.dy  -  au^b^ 81^  ydy  -  Uj„ifb/'8i  dy 
o'  o  1  o'  o  1 

Hg  =  ab^ /’'''82bydy  +  r?b/'82bdy  -  au^b^  /  '  Sj ^ydy  -  u„r?b 8u,dy 


M  .H  ,  9  -H 


H4  =  ah^  f  S2Kydy  +  rfbf  S2Kdy  -  au^b  f  8ij^ydy  -  “111^?^/  S^j^dy 
o 

‘/2 


Ha 


-2  r  c  ,,J,,  ,.0!k  r  Q  A,r  J" 


ah^f  S2^ydy  -  r?b  /  82xdy 


„a+l 


II  =  «b2/'83pydy  +  r?b  /  83pdy  +  ^  (pjUj) 


U  =  ah^  f  83^, ydy  +  r?b  /  83^,dy 


13  =  ah^  /^Sgj^ydy  +  rfb  /  83j,dy 

O  O 

14  =  ab^  f  S3j^ydy  +  r'Jb  /  83j,dy 

O  o 

15  =  -  a.b2/^83.^ydy  -  rfb/^83^dy 

o  O 
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The  wake  regime  coefficients  (F,  G,  H,  I,  and  J),  where  g  =  Uc ,  rj,  b,  K,  Cc  and  x  are 


,a+l 


Fi  =  ydy  +  r?b/'si,  ^7  + 

o  c  O  ° 

F2  =  ah^  I  Si^.ydy  +  rfb  /  S^^.dy 

0  1  0  1 


Fg  =  ab^/  Sgj^ydy  +  r?b  /  S^^dy 

0  0 

F4  =  ab^/  Sjj^ydy  +  r?b  /  Sjj^dy 


Fg  =  ah^  f  Sgc  ydy  +  r?b  /  dy  + 

o  c  o  ^ 


.a+1 


{a+  1)  dC^ 


(Pr^r^ 


^6  ^ix'^y  -  ps^s^'s 


9  1  1  r9  o 

Gl  =  «b  /  Sg^^ydy  +  rfhf  Sg^^dy  ^(Pc'^J 

Go  =  a  b^  f  So_  ydv  +  r?b  f  Sr.  *^7 
z,  •'  1  J  zr. 


Gg  =  ab^f  Sggydy  +  rfb  /  Sggdy 


G4  =  ab^  f  S2j.ydy  +  r?b  /  Sgj^dy 


9  1  1  d  2 

.2  r  c  -_j..  .  _ai.  r  c  J  1  ^  (p  u^) 

'■r  r  c  ^ 


Gg  =  ab^f  Sgc  ydy  +  rfb  f  Sg^  dy  +  ^ 


(a  +1)  dC 


"6  =  /  ^2x7^7  -  r?b/  Sg^dy  -  p^u^v^r 


Hi  =  ab^/^Sg  ydy  +  rfb 82^  dy  -  u^ab^/^'Si^  ydy  -  u^rfb/sj^  dy 


,a+l  ,a+l 

i  ^  /o  2n  i  (9 

(Pc'^c>  -  '^mT- 


(a  +  1)  c 


^(a+  1)  (9u^ 


Hg  =  ab^/'Sg^.ydy  +  rfb/'^Sg^  dy  -  b^  / '  S  j ^  ydy  -  u^r?b Sj^  dy 
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H3  =  ab^/  S2bydy  +  rfb/  S2i,dy  -  u^a  b^  /  S^^ydy  -  u^rjb  /  S^^dy 


H4  =  a  b^ S2j,ydy  +  r'?b Sj^dy  -  V  S^^ydy  -  u^r?  b S^j^dy 


H5  =  ab^/'Sjc  ydy  +  rfb/  Sjc  dy  -  u^ab^  f  Sj^  ydy  -  u^rfb  /  dy 


r«+^  .  r“+^  . 

1 _  o  /  2s  1  o  , 

(a+1)  dC^  ''m(„  +  1) 


Hg  =  -  ah^  S2^ydy  -  rfb/'  S2^dy  +  u^a  b^ /'  S^^ydy 


+  S^^dy  +  vr 


r«+l 

d 


II  =  /  S3  ydy  +  rfbf  ^  (Pc'^cCJ 


I2  =  ab^  f  Sg^.ydy  +  rf  b  /  S3^,dy 


13  =  ab^f  Sgj^ydy  +  rf b  /  Sgj^dy 

o  o 

14  =  ab2  /^Sgj^ydy  +  rfb  f  Sgj^dy 

o  O 

1  1  '■?+^  n 

h  =  ^30  y'ly  +  ^3c  +(^^) 


Ig  =  -  a 


Sg^ydy  -  rfb/  Sg^dy 


Jl  =  al>2/'S3^  ydy  +  rfb  /  '  Sg^  dy  -  C^ab^/'Si^  ydy  -  C^rfb/’'“Ss  dy 
0^  O^  0^  oc 

rCt+l 

^  (Pc’^c^c^  -  •^mT^rTTY 


(a+1)  -  -m(„+l) 


J2  =  ab^/  Sg^  ydy  +  r?b  /  Sg^  dy  -  C^a  b^  /  '  Sj^  ydy  -  C^^rPb Sj^  dy 

O  1  n  1  O  1  n  1 
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Sa^ydy  +  rfb S^^dy  -  c^ab^/'  Sjj^ydy  -  C^rf b  /  Sjj^dy 

O  O  O  O  — ~ — 

=  Sg^ydy  +  rfb/'"''  Sgj^dy  -  C^a  b^ /'  Sjj^ydy  -  C^r?b/'si„dy 


J4 

J5  =  ydy  +  r?b7^^S3c  dy  -  C^ab^/"Si(.  ydy  -  C^rfb /"  dy 


,  ■  .‘/2 


[.a+1  rfl+l 

W«„cj  -  c„niT-r,  4- (p„"„) 


(o  +  1)  ''-c“c“c 


'm(a  +  1)  dC^ 


Jg  =  -  ah^ S^jdy  ~  rfb/'Sg^dy  +  C^a  b^ /'  Sg^ydy  +  C^rfb /' Sg^dy 

000  O 


+  Q  r 

r 


In  the  fully  developed  regime  define 


y  =  r/b 

r  =  b 
s 

r“dr  =  b(by)®  dy 

The  fully  developed  regime  coefficients  (F,  G,  H,  I,  and  J),  where  g  =  p, Uc,  b,  K,  Q  and  x  are: 

Fl  =  bb«/'Sipy«dy 

o 

Fg  =  bb«  /  y«dy 

o  C 

Fg  =  bb«/  Si,^y«dy 
0 

F4  =  bb«/  Sjj^y«dy 

0 

F5  =  bb«/'sic  y“dy 

o  C 

Fg  =  bb«/  Sj^y“dy  _  PgVgr« 
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Gi  =  bb«/  S2py«dy 


^a+l 


(a  +  1) 


^2  =  bb«  /  S2^  y«dy 


0  c 


Ga  =  /  S2by“dy 

0 

G4  =  bb“/  S2j(.y“dy 


Gg  -  bb“  f  S2Q  y'^dy 


’2C 

o  c 


Gfi  =  -  bb«/  S2  v«dy  - 


p  V  u  r 
^  s  s  s  s 


.a+1 


Hi  =  bb« S2py«dy  -  u^bb«/^'Sipy«dy  + 


H2  =  bb«/'S2„  y“dy  -  u^bb“/'Si^  y«dy 


Hg  =  bb«/^S2i,y“dy  -  u,„bb«/^Sigy«dy 
o  O 

H4  =  bb«/'S2Ky“dy  -  %bb«/''''  Sn^y“dy 


'72 


Hr  =  bb«  /  S2C  y“dy  -  u^bb«  /  Sgc  y«dy 


H, 


bb«/'S2j“dy  +  u^bb“/  Sj^dy  + 


11  =  bb“/  S3  y«dy 

0 

12  =  bb«/'s3^  y“dy 


I3  =  bb«/  Sgi^y^dy 


[4  =  bb“/  Sgj^y^dy 
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I5  =  bb«/'s3c  v% 

O  C 

le  =  -  S3^y“dy 

O 

Jl  =  bb^/'Sg  y^dy  -  C^bb“ Si  y^'dy 
o  0 

J2  =  ^3^  y«dy  -  c^bb«/*'''  Sj^  y«dy 

o  ^  o  ^ 

J3  =  bb®  /  '  Sgi^y^dy  -  C^bb« /' Sji^y^dy 
o  o 

J4  =  tb«/\S3^y«dy  -  c^bb“ Sij,y«dy 

o  o 

J5  =  y“dy  -  C„bb«/“Si(i  y“dy 

o  ^  o  ^ 

Je  =  -  S3^y«dy  +  C^bb«/^'Si^y«dy  +  q^r« 

o  o 
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NOMENCLATURE 

b  Shear  layer  thickness 

Cp  Specific  heat  at  constant  pressure 

C  Element  fraction  of  bleed  gas 

F,G,H,I,J  Coefficients  in  system  of  ordinary  differential  equations 
g  Generalized  variable  (Appendix  A) 

h  Base  height  in  planar  flow 

K  Parameter  in  species  profile  equation 

M  Mach  number 

n  Boundary-layer  profile  exponent 

N  Generalized  coefficient  in  system  of  ordinary  differential  equations 

p  Static  pressure 

pb  Base  pressure 

Po  Total  pressure 

Po'  Pitot  pressure 

P  Pressure  gradient  parameter 

q  Turbulent  species  flux 

r  Lateral  coordinate 

rb  Base  radius 

q  Lateral  distance  from  centerline  to  inner 

rs  Lateral  distance  from  centerline  to  outer 

R  Gas  constant 

RSP  Rear  stagnation  point 


edge  of  shear  layer 
edge  of  shear  layer 
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Ree 

Momentum  thickness  Reynolds  number 

Rx 

Turbulent  Reynolds  number 

S 

Generalized  integrand  (Appendix  A) 

To 

Total  temperature 

u 

Axial  velocity  component 

V 

Lateral  velocity  component 

X 

Axial  coordinate 

Ax 

Integration  step  size 

y 

Distance  from  wall  in  initial  boundary  layer;  dimensionless  mixing  zone 
coordinate  (Appendix  A) 

!  a 

Geometric  parameter  (one  for  axisymmetric  flow,  zero  for  planar  flow) 

1 

Ratio  of  specific  heats 

5 

Initial  boundary-layer  thickness 

Distance  at  beginning  of  recompression  between  centerline  and  streamline 
which  orginates  at  outer  edge  of  initial  boundary  layer 

6 

Momentum  thickness  of  initial  boundary  layer 

P 

Density 

T 

Turbulent  shear  stress 

SUBSCRIPTS 

1 

Conditions  before  base  expansion 

a 

Outer  free  stream 

c 

Centerline 

d 

Conditions  on  dividing  streamline 

j 

Conditions  in  base  bleed  flow 
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m  Conditions  at  middle  of  shear  layer 

p  Conditions  at  beginning  of  recompression 

r  Conditions  at  rear  stagnation  point 
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